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The role of mixed state entanglement in liquid-state nuclear magnetic resonance (NMR) quantum 
computation is not yet well-understood. In particular, despite the success of quantum information 
processing with NMR, recent work has shown that quantum states used in most of those experiments 
were not entangled. This is because these states, derived by unitary transforms from the thermal 
equilibrium state, were too close to the maximally mixed state. We are thus motivated to determine 
whether a given NMR state is entanglable - that is, does there exist a unitary transform that 
entangles the state? The boundary between entanglable and nonentanglable thermal states is a 
function of the spin system size N and its temperature T. We provide new bounds on the location 
of this boundary using analytical and numerical methods; our tightest bound scales as N ~ T, 
giving a lower bound requiring at least N ~ 22, 000 proton spins to realize an entanglable thermal 
state at typical laboratory NMR magnetic fields. These bounds are tighter than known bounds on 
the entanglability of effective pure states. 

PACS numbers: 03.67.Mn, 03.67.Lx, 76.60.-k 



I. INTRODUCTION 

Entanglement is hidden information in the form of non- 
local correlations between two or more quantum degrees 
of freedom. Recent work supports the view that entan- 
glement is a quantifiable physical resource. Entangled 
states are useful for communicating quantum informa- 
tion 0,0) HI- I* also appears that entangled states are re- 
quired to exploit the exponential speedup of many pure- 
state quantum algorithms 0,0, including Shor's algo- 
rithm BULB- 

Quantum computers based on liquid-state nuclear 
magnetic resonance (NMR) have successfully imple- 
mented these same quantum algorithms 0, 0, 0> 
Il3j . However, Braunstein et al. have shown that the ex- 
periments use separable states The initial NMR 
states are not pure, but rather mixed states derived from 
a thermal ensemble. Due to the small number of qubits 
and the weak initial polarization of the ensemble, these 
effective pure states are too random to possess any en- 
tanglement. 

This startling result raises doubts about the validity of 
NMR quantum computation. How can NMR techniques 
demonstrate quantum algorithms without entanglement? 
Are entangled states necessary resources for quantum 
computation? These questions are yet to be conclusively 
answered, although it has been pointed out that 1) the 
essential role of entanglement in pure-state algorithms 
does not necessarily extend to the mixed states used in 



'Current address: Department of Physics, Yale University, New 
Haven, CT 06520. 



NMR 0, and that 2) there exist efficient pure-state 
quantum algorithms without entanglement [l6j. In fact, 
a model of computation using one qubit purified from a 
mixed state has been shown to be capable of performing 
non-classical algorithms Furthermore, attempts to 

construct classical models of NMR quantum computation 
have thus far suffered from exponential scaling deficien- 
cies 0, 0] . Some speculate that entangling operations 
— rather than entangled states — are the source of ex- 
ponential speedup in quantum computers pol |. 

Anwar et al. recently implemented the Deutsch-Jozsa 
and Grover algorithms on a nearly pure two-qubit NMR 
state 0, , proving that it is possible to perform NMR 
quantum computation with entanglement. These results 
are impressive, but it will be challenging to scale the sys- 
tem since each qubit must be selectively excited. More- 
over, as the authors state |23| . the experiment cannot an- 
swer fundamental questions regarding the role of mixed 
state entanglement in quantum computation. 

We attempt to address these questions through a theo- 
retical study of entanglement in thermal states, the nat- 
ural mixed states in liquid-state NMR systems. Due to 
reasons we explain, the thermal state might be more eas- 
ily entangled than the effective pure state used in current 
NMR quantum computation experiments |53| . 

The spins in a liquid-state NMR system are largely 
uncoupled, although the tiny coupling allows a universal 
set of unitary operations to be performed on the spins. 
We thus approximate the initial state of the NMR sys- 
tem with a fully uncoupled Hamiltonian. For ease in our 
analysis, we also assume that all the spins are identical. 
The thermal state corresponding to such a system is then 
completely specified by the number of qubits (spins) N 
and the dimensionless quantity a, which is a measure of 
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the state's polarization. The polarization is the bias of 
the spins towards parallel alignment with the magnetic 
field; it is also approximately invesely proportional to the 
equilibrium temperature T. Finally, we assume that an 
arbitrary unitary transform can be efficiently performed 
on each N qubit system in the ensemble. 

Here we present bounds on the boundary in N-a pa- 
rameter space between thermal states that can and can- 
not be entangled. A lower/upper bound on the bound- 
ary demarcates parameter space where thermal states 
can/cannot be entangled. We obtain a lower bound on 
thermal states that is tighter than previous calculations 
for effective pure states 1^ ■ From a practical stand- 
point, a lower bound on the parameter space boundary is 
equivalent to an upper bound on the minimum number of 
qubits and polarization required to create an entangled 
NMR state in the laboratory. 

We begin in Sec. [H] with a brief overview of previous 
work on NMR state entanglement. Then we assess differ- 
ent methods for entangling NMR states in Sec. II I II and 
describe our approach to calculating and bounding the 
entanglement of transformed thermal states in Sec. IIVI 
Next, we present the results of this paper. In Sec. [3 we 
use the negative partial transpose test [2^, |2(j to calcu- 
late bounds on the entanglement of thermal states trans- 
formed by a specific Bell unitary. In Sec. I VII we develop a 
general method to calculate bounds on the entanglement 
of generic Bell-transformed thermal states, based on a 
classification scheme by Diir and Cirac |24|. Combining 
the method with majorization theory could potentially 
lead to tighter bounds. This possibility is discussed in 
Sec. lVIll Finally, we summarize our results and conclude 
in Sec. IVTTT1 

II. REVIEW OF BOUNDS ON NMR STATE 
ENTANGLEMENT 

We first briefly review mixed state entanglement. A 
mixed state p shared by party A and party B is separable 
if and only if it can be expressed as 

i 

where pf and pf are density matrices contained in the 
Hilbert spaces of parties A and B respectively and the 
weights Pi are probabilities. The bipartite state p is en- 
tangled if it cannot be written in this form. 

In general, a mixed state p shared by K parties is sep- 
arable if and only if it can be expressed as 

P = 5>®"?- (2) 

i j=l 

Each density matrix p\ is in the Hilbert space of the jth 
party. The A-partite state p is entangled if it cannot 
be written in the form of Eq. If A' is equal to the 



number of degrees of freedom (e.g. the number of qubits) 
and p satisfies Eq. Q, then we call p fully separable. 

We desire a way to determine whether a state is entan- 
gled. The above definitions for entanglement are a pos- 
sible starting point, but they require us to examine an 
infinite number of decompositions. For bipartite mixed 
states, there exists a simple, computable criterion for en- 
tanglement 25, 26]: 

If the partial transpose of a bipartite state p 
has a negative eigenvalue, p is entangled. 

Note that this criterion gives only a necessary condition 
for entanglement except in the case of two qubits, where 
it becomes both necessary and sufficient. 

The partial transpose operation is defined as follows. 
Any bipartite density matrix p may be expressed as 

p= C a,a',b,b'\a)(a'\(g>\b)(b'\. (3) 

a,a' ,b,b' 

The basis states \a), \a') and \b), \b') are in the Hilbert 
spaces of parties A and B respectively, and the numbers 
C a .a>.b,b> are complex and constrained to give p unit trace. 
The partial transpose of p with respect to party A is 
denoted p TA and given by 

P TA = J2 C *-*'-b,b>W){a\®\b){b'\. (4) 

a,a f ,b,b f 

If p TA has at least one negative eigenvalue, we say that p 
has negative partial transpose (NPT). Otherwise, we say 
that p has positive partial transpose (PPT). 

We now turn to the specific case of entanglement in 
NMR quantum computation. Liquid-state NMR quan- 
tum computing experiments are performed on samples 
containing an ensemble of 10 18 — 10 19 identical molecules 
under a strong, static magnetic field (~12 T). Each 
molecule possesses N spin- 1/2 nuclei, which represent 
the qubits. The large number of molecules is needed be- 
cause the signal from a single nuclear spin is exceedingly 
weak. The dynamics of the NMR system is dominated 
by the Zeeman interaction between the nuclear spins and 
the magnetic field. 

The natural quantum state in NMR is the thermal 
state 

e -H/kT 

Pth = — -g — (5) 

where Ti. is the Hamiltonian, k is the Boltzmann constant, 
T is the absolute temperature, and Z is the partition 
function. 

The thermal state is highly mixed in liquid-state NMR 
systems. Yet most quantum algorithms require pure 
states. To circumvent this problem, experimenters first 
create initial spin states of form 

p ee =(l-e)M d + e\0){0\. (6) 
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Here the dimension of p e g is d = 2 N , = Ig/d is the 
maximally mixed state (Id is the d-dimensional identity 
matrix), and e characterizes the fraction of ground state 
1 0) . These states are called effective pure states or pseu- 
dopure states [27ll28j because 1) a unitary transform only 
evolves the excess ground state population and 2) the 
term proportional to Md is unobservable in conventional 
NMR experiments. 

Braunstein et al. established that the states used in all 
NMR quantum computing experiments up to now have 
never been entangled [bH ]. They studied near maximally 
mixed states of the form 



p e = (1 - e)M d + ep' 



(7) 



where p' is an arbitrary density matrix and e character- 
izes the fraction of the mixture in p' 54] . Effective pure 
states clearly fall into this class of states. 

The authors established bounds on the entanglement 
of p e for a system of N spin- 1/2 particles: 



e < 



e > 



1 



1 + 2 2iV-l 
1 

1 + 2 N / 2 " 



=>■ p e is always separable (8) 
p e can be nonseparable . (9) 



The first expression may be interpreted as a lower bound 
on the size of the separable neighborhood around M d - 
When e is sufficiently small, p e is separable. The sec- 
ond expression may be interpreted as an upper bound on 
the separable neighborhood. When e is sufficiently large, 
there exists some p' such that p e is nonseparable. 

These bounds are applicable to liquid-state NMR be- 
cause quantum computations can only transform a state 
of form p € to another state of form p t . Liquid-state NMR 
systems are approximately closed and thus evolve unitar- 
ily. Assuming the largest practical initial polarization, e 
may be as large as 3 x 10 -5 for N=2. Since e scales as 
N/2 N for low polarization, Eq. JSJ shows that the states 
used in NMR quantum computation are separable when 
N < 12. Therefore, current NMR experiments, for which 
N < 7, do not use enta ngle d states. 

Gurvits and Barnum [29j recently tightened the lower 
bound to 



e < 



2(6 W / 2 



(10) 



This result shows that the states used in NMR quantum 
computation are separable when N < 32. 

There is also a tighter upper bound due to Diir and 
Cirac, but we defer discussion of this bound to Sec. IVII 

In what follows, we will be interested in NMR param- 
eter space described by N and a (defined in Sec. IHIj l. 
We thus re-express the above bounds in terms of these 
parameters. 

It can be shown that e is given by 



/z 



2 N -1 



Na 



(11) 



with the approximation being valid in the limit a <C 1. 

Using Eq. (|11|) . the Braunstein et al. bound on en- 
tanglable p e g is 



a>~]n(V2-l), 



(12) 



and the Gurvits-Barnum bound on nonentanglable p e R is 



a < 



3(2 N - 1) 
2^(6^/2) 



(13) 



In this paper, we call a state entanglable if there exists 
a unitary that transforms the state into an entangled 
one. This terminology allows us to discuss entanglement 
without referring to the specific unitary that entangled 
the original state. 



III. HOW TO ENTANGLE AN NMR STATE? 

We wish to address the question of whether it is pos- 
sible to entangle an NMR thermal state. In this section, 
we evaluate several methods for achieving entanglement. 
The discussion here motivates our approach in Sec. IIVI 

A. Initial NMR state 

The first consideration is the choice of initial state. 
Current NMR quantum computation begins with effec- 
tive pure states, but we simply want to obtain an entan- 
gled state by any means. The thermal state is a better 
initial state for reasons we now explain. 

One possible decomposition of an effective pure state 
is the following mixture of transformed thermal states: 



d-l 



Pcft - 



^PjUjPtbtf . 



(14) 



3=1 



Here d is the dimension of p e s, pj = l/(d — 1), and Uj 
are cyclic permutation matrices. 

Since the numbers pj are probabilities, p c ff satisfies the 
mathematical definition of a convex combination. Now 
entanglement is generally a convex function and there- 
fore is reduced under convex combination. Hence, an 
individual transformed thermal state UjPthUj may pos- 
sess more entanglement than an effective pure state of 
the same dimension. 

These considerations motivate us to focus on entan- 
gling NMR thermal states in this paper. 

For convenience in the analysis, we assume that all 
spins in the NMR molecules have the same Zeeman en- 
ergy splitting and consider an isotropic Hamiltonian 



(15) 
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expressed in the computational (spin) basis. Here v = 
7_B is the frequency corresponding to a spin with gyro- 
magnetic ratio 7 in a magnetic field of strength B, Zi 
is the Pauli Z operator acting on the ith spin, and N is 
the number of spin-1/2 particles in each molecule. We 
neglect the J-coupling interaction between neighboring 
spins as it is 10 -6 times smaller than the Zeeman energy. 

The corresponding thermal state is given by Eq. (JSJ) 
and the Hamiltonian above. Notice that p t h is diagonal 
in the computational basis, invariant under the exchange 
of any two spins, and fully separable. 

We define a, a measure of the thermal state's polar- 
ization, to be 



Notice that when a -C 1, this quantity is approximately 
the difference between the fraction of spins aligned and 
the fraction of spins anti-aligned with the magnetic field, 
i.e. the polarization. Since NMR quantum computing ex- 
periments are performed in the a -C 1 regime, we hence- 
forth refer to a as simply the polarization in this paper. 

Putting together Eqs. JSJ, l|15l) . and (|16[l . the matrix 
elements of the thermal state in the computational basis 
are 

(i\p^\j) = |^- 2 -«]« (17) 

where the Hamming weight w(i) is the number of Is in 
the binary expression for i and < i,j < 2 N — 1. In the 
limit a — > 0, we see that pth - * I- 

B. Methods to facilitate entanglement 

An experiment to entangle an NMR thermal state con- 
sists of two parts: 

1. Prepare the initial thermal state pth with number 
of qubits N and polarization a. 

2. Apply a unitary operation to entangle the thermal 
state: p th >-> Up^W. 

In view of the above procedure, there are several ways 
we can facilitate the creation of entangled thermal states: 

• Enhance initial polarization 

• Increase the number of qubits (spins) per molecule 

• Perform algorithmic cooling 

• Find optimally entangling unitary operations. 

We expect that increased polarization will yield a ther- 
mal state that is more easily entangled. Doing so moves 
the thermal state away from the maximally mixed state. 
Increasing the number of qubits should have the same 
effect, but the intuition behind this conjecture is weaker. 



As we have seen, the bounds on the nonentanglability 
of p e g relax as N increases. Moreover, the volume of 
separable states appears to decrease with Hilbert space 
dimension. A lower bound for the volume of separa- 
ble states [3(J has been found, which decreases exponen- 
tially with N. Numerical evidence |3j supports the same 
trend. 

We now examine the viability of each method. 



1. Enhancing initial polarization 

We discuss three experimental approaches for increas- 
ing the initial polarization: 1) lowering the temperature, 
2) increasing magnetic field strength, and 3) chemical and 
optical methods. As a reference, a state-of-the-art liquid- 
state NMR experiment using seven proton nuclear spin 
molecules at room temperature has parameters N = 7, 
q = 4x 1CT 5 , B = 12 T, and T = 300 K. A seven-qubit 
effective pure state must have a > 3 x 10~ 2 to be outside 
of the nonentanglable regime described by Eq. I|13l) . 

Lowering the temperature of the experiment is imprac- 
tical because any appreciable improvement in polariza- 
tion demands such a dramatic drop in temperature that 
the sample becomes solid. Unless the spin system is di- 
lute, it acquires large dipolar couplings that smear out 
the spectrum and make selective excitation of spins diffi- 
cult. A dilute spin system could be used, but we ideally 
desire large N. There are several proposals for solid- 
state NMR quantum computers that address these prob- 
lems |32l l33| , but none have been experimentally realized 
thus far. 

Using a stronger magnetic field as a way to increase 
polarization is also limited. Because present NMR mag- 
nets are based on superconducting coils that are current- 
limited, the largest attainable field is 21 T. Higher mag- 
netic fields can be achieved with DC Bitter resistive mag- 
nets (33-45 T) and capacitively-driven magnets (50-60 
T), but the generated fields are too spatially inhomoge- 
neous to be used in liquid-state NMR. The magnets also 
consume massive power and require extensive cooling. 

Inducing higher polarization by chemical and opti- 
cal means is a much more promising approach and has 
already been experimentally demonstrated. Optically 
pumped xenon has been shown to give a tenfold increase 
in a for a two qubit molecule |34[ ■ Another factor of three 
or four may be gained if the xenon is isotopically pure. 
Further increases in a are also possible if higher polar- 
izations of xenon are used. The experiment of Ref. |34| 
used 1% polarized xenon, but Xe polarizations as high 
as 67% have been obtained j3f|. Methods involving para 
hydrogen (two protons in the singlet state) yield even 
larger polarization. By reacting para hydrogen with an- 
other molecule, Anwar et al. [23| have recently created 
two-qubit effective pure states with e = 0.916 ± 0.019. 
These initial states are entangled according to Eq. ijjjj . 
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2. Increasing the number of qubits 

Entanglement can also be achieved by using more 
qubits, but this approach is limited by substantial fun- 
damental problems |36j . Since individual spins are ad- 
dressed by selective excitation, the resonance frequency 
of each spin must be well-separated from the others. This 
requirement may be satisfied by using nuclear spins of 
different chemical species. There are five distinct spin- 
1/2 nuclei commonly used for liquid-state NMR: 1 H, 13 C, 
15 N, 19 F, and 31 P. Spins of the same type may also 
be used, but the largest number of homogeneous spins 
demonstrated in an experiment thus far 37] is N = 6. 
Thus, for conventional liquid-state NMR quantum com- 
putation, the requirement of selective excitation restricts 
experiments to N — 30. It may be possible to im- 
prove upon this limit by using polymers |3Sj or higher- 
order spins. Classical logic operations have been experi- 
mentally demonstrated in spin-3/2 |39j and spin-7/2 nu- 
clei E3. 



3. Algorithmic cooling 

Algorithmic cooling is a compression technique that 
concentrates the polarization of many qubits into a 
smaller number of qubits. A well-known example is the 
Schulman-Vazirani procedure |4l| . which starts with Nq 
qubits in thermal equilibrium at polarization ao and ex- 
tracts from these initial qubits a string of order (ao) 2 Ao 



qubits each with po > 1 



2 N 



-10 



where po is the proba- 



bility of measuring |0). This method produces very pure 
qubits but at a high cost for the weak polarizations in 
current NMR experiments. For «o = 4 x 10~ 5 , we need 
approximately 6.3 x 10 8 initial thermal qubits to obtain 
just one pure qubit. The enormous number of initial 
qubits is needed because the Schulman-Vazirani proce- 
dure preserves the entropy of the system. If entropy 
conservation is circumvented, better results are possible. 
Another scheme proposes to use an external set of spins 
that act as a heat bath ji^. This method can create 
effective pure states with initial to final qubit ratio of 50. 

If a can be boosted with some of the techniques pre- 
viously mentioned, the Schulman-Vazirani procedure be- 
comes more realistic. For example, beginning with No = 
30 and ao = 0.2, we can obtain one nearly pure qubit. 



Qubit 1 
Qubit 2 
Qubit 3 
Qubit N 





H 




— (— 
















V 


J 












• 
• 


r 




r 


>- 









Ufa 



FIG. 1: Quantum circuit for [/cH^fan unitary operation. The 
I7f an unitary acts like a fanout gate on computational states. 
The C/cH unitary transforms computational states into Bell 
states. 



IV. APPROACH 

We choose to investigate the approach of applying a 
unitary transform to entangle the thermal state. More 
specifically, we are interested in the question of whether 
a given iV-qubit thermal state with initial polarization 
a is entanglable, that is, does there exist a unitary that 
entangles the state? 

Studying this problem will be fruitful for two reasons. 
First, we can address fundamental questions regarding 
NMR state entanglement by examining the scaling be- 
havior of entanglable regions in N-a parameter space. 
Second, our study provides insight as to how experiments 
may be designed to realize entangled NMR states. 

To establish whether a thermal state is entanglable, we 
must in principle search over the space of all unitary oper- 
ations. Rather than attempt this intractable procedure, 
we focus our attention on the family of transformations 
C/b that maps the computational basis to the Bell state 
basis. We choose Ub as a starting point because it gener- 
ates maximal entanglement when applied to the ground 
state. The thermal state has the greatest population in 
the ground state, so Uh may be effective in entangling p t h- 
As we shall see, the Bell unitaries also possess symmetries 
that make them particularly amenable to analysis. 

We first define a "standard" Bell unitary: the con- 
trolled NOT-Hadamard transform U c k depicted in the 
right half of Fig. ^ A Hadamard gate is applied to the 
first qubit followed by a set of controlled-NOT (CNOT) 
gates on the rest of the qubits. The other unitaries in 
the Bell transformation family are permutations of U c n: 



4- Entangling unitary operations 

Entangling NMR states via unitary operations has 
been largely unexplored. The thermal state is separa- 
ble, but we can apply a unitary transform to move pth 
into an entangled region in Hilbert space. It is currently 
unknown what unitaries optimally entangle mixed states 
except for the case of two qubits Q • One group has pro- 
posed quantitative measures for the power of entangling 
operators [44| . 



U h (P) = Pt/cHPt (18) 

with P being a matrix that permutes the computational 
state to Bell state mapping. 

We will later study another Bell unitary U c nUf an , 
shown in Fig. ^ The permutation matrix Uf an is also 
known as the fanout gate because it bit flips the latter 
N — 1 qubits if the first qubit is in state |1). 

For notational convenience, we also define the following 
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transformed states: 



A. Negativity formula for {1,N — 1} bipartite split 



p h = U h p t hUl 
PcH = U c npthUl H 

Pel = UdiUfanPthU^U^ 



(19) 
(20) 
(21) 



Our task is now to determine if a given Bell- 
transformed thermal state pb is entangled. One method 
is to apply the negative partial transpose test. This ap- 
proach is studied in Sec. El 

If all the eigenvalues of partial transposed state can be 
calculated, we may quantify the amount of entanglement 
by the negativity entanglement measure [3ll |45| . The 
negativity of a state p is defined as 



EM 



\P TA h 



1 



(22) 



Here \\p TA ||i denotes the trace norm of p TA , which is de- 
fined as tr(p TA<! p TA ). For diagonalizable operators, the 
trace norm reduces to the sum of the absolute eigenvalues 
J2i \W- Negativity has the range < Ej^{p) < 1/2, tak- 
ing on the smallest value for separable p and the largest 
value for maximally entangled pure states. Compared to 
other known entanglement measures 0, S3 > it has the 
advantage of being relatively simple to compute for arbi- 
trarily large N. 

The partial transpose test may give different results de- 
pending on how the qubits are distributed between par- 
ties A and B. The thermal state is totally symmetric, 
and the unitary transforms we examine in this paper are 
symmetric among the latter N — 1 spins. Therefore, it 
is sufficient to consider partitions where the first q spins 
are assigned to party A and the remaining N — q spins 
to party B. We write such a bipartite split as {q, N — q}. 

A more general method for determining the en- 
tanglability of thermal states is the following. Suppose 
we have a transformed thermal state whose entanglement 
is difficult to determine, but we also have a state p' whose 
entanglement is easy to classify. If we can find a quantum 
operation that converts the transformed thermal state 
into p' without generating any new entanglement, then 
the entanglement of p 1 bounds the entanglement of the 
transformed thermal state. This approach is studied in 
Sec. EH 



V. NPT BOUNDS ON THE ENTANGLEMENT 
OF CNOT-HADAMARD-TRANSFORMED 
THERMAL STATES 

We first derive a formula for the negativity of the 
CNOT-Hadamard-transformed thermal states under the 
{1, TV— 1} bipartite split. Then we find an NPT bound on 
entangled CNOT-Hadamard-transformed thermal states 
under the {N/2, N/2} bipartite split. This bound is em- 
pirically derived, but it will be confirmed by the analysis 
of Sec.EU 



The quantum circuit for U c h (see Fig.QJ suggests that 
the first qubit is special. The thermal state may be 
rewritten to separate the behavior of the first qubit from 
the others: 



Pth = exp 



exp ( — Ziaj <g> exp 




(23) 



where we have used Eqs. (jSJ, 1|15(1 . and (IB) - Operators 
with a subscript 1 are taken to be acting on the first 
qubit. 

We then insert a complete set of angular momentum 
states \j, m) (j, m\ to represent the collective spin state of 
the latter N — 1 qubits. The symbol j denotes the total 
spin angular momentum whereas m = — j, + — 
1, j denotes the total azimuthal spin angular momentum. 
As a result, the above equation becomes 



E-2ma 

3,m 



j,m)(j,m\ 



(24) 



Because (^2^ =2 \j, m) = m\j, m), each angular mo- 

mentum term is multiplied by e _2ma . 

Now we examine the effect of U c n on the thermal state. 
It is easier to analyze the problem if we break the unitary 
down into the application of a Hadamard gate on the first 
qubit (Hi) followed by a collective CNOT on the other 
N - 1 qubits (C/cnot), i-e. J7 c h = U cnot Hi. 

Applying a Hadamard operation to the first qubit of 
pth yields the state 



Pth = 



H lPth Hl 
1 



E-2ma 



\j,m)(j,m\ (25) 



since Hi = H\ and H x ZiH x = X x . Expanding e Xia , 
this equation can be rewritten as 



Pth = -= ^cosha)/j - ( siiiho ) .Yj 



-2rna 



\j,m)(j,m\ 



= i[cosha(|0)(0| + |l)(l| 
+ sinha(|0)(l| + |l)(0|)j <g)^\ 



(26) 

- 2ma \j,m){j,m\. 



In the last line, we have expanded I± and X x in the com- 
putational basis of the first qubit. 

To obtain the full transformed state p c u, we apply 
the collective CNOT. Since U cnot \0) |j, m) = \0)\j,m) and 
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FIG. 2: (color) Negativity (a) and logarithm of negativity (b) of CNOT-Hadamard-transformed thermal states under {1, TV — 1} 
bipartite split in TV — a parameter space. The logarithm of negativity is calculated as log 10 \Ej^(p c n) + 10~ 12 ] . Transformed 
thermal states p C H are nonseparable in the parameter space left of the bound marked "p C H nonseparable" (Eq. 1)31(1 ). Effective 
pure states are entanglable in parameter space left of the Braunstein el al. bound (Eq. 1121 ) and are nonentanglable in the 
parameter space beneath the Gurvits-Barnum bound (Eq. 113H ). 



C^cuot 1 1) I J j Tn) = \l)\j,—m), we see that p C H is composed 
of two-dimensional subspaces spanned by \0)\j, to) and 
\l)\j,-m): 



-2ma 



cosh a sinh a 
sinh a cosh a 



(27) 



Next, we calculate the partial transposed state p^ 
with respect to the {1, N— 1} split. The partial transpose 
solely mixes the states in the subspace of \0)\j, ±to) and 
l 1 )b')±"^)» giving 



-2ma , 



(28) 



Each term of the above sum has eigenvalues 

e~ 2ma cosh a ± e 2ma sinh a 
A±(mJ = 



(29) 



The largest possible m-value is (N — l)/2, so the min- 
imum eigenvalue of is 



An 



e -(W-l)a coshQ ,_ e (W-l)a sinh( 



(30) 



When A m in is negative, p C H is entangled according to the 
negative partial transpose test. Since A m i n is monotonic 



in N and a, the equality A m i n = provides a bound on 
p C H being nonseparable: 



e 2(A '- 1 ) Q tanha> 1. 



(31) 



As Eq. specifies all the eigenvalues of , we also 
calculate the negativity by summing A± over all j and 
to- values: 



-EV(Pch) 



-1 



Ej, m (|A+(m)| + |A_(m)| 



1 1 



2j D 



2 2 \ Jmax -J J j max + j + 1 



2j + l 



(|A+(m)| + |A_(m)| 



(32) 



m=—j 



Here the maximum j- value is j ma x = (N — l)/2. The 
factor in front of the sum over to is the multiplicity of j, 
which is equal to the number of states with m = j minus 
the number of states with m = j + 1. For N — 1 even, 
the range of j is 0, 1, . . . j ma x, whereas for N — 1 odd, the 
range of j is 1/2, 3/2, . . . j max . 

Fig. 2(a) shows the negativity plotted up to -/V=500. 
Note that the horizontal axis is plotted in units of 
log 10 (a _1 ) such that smaller polarization (weaker mag- 
netic fields and higher temperatures) lies towards the 
right. As expected, the minimum a required for an en- 
tangled thermal state decreases with N. However, the 
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slope of negativity around this a becomes increasingly 
shallow. This behavior is seen most clearly in Fig. |2 (b)| 
where we have plotted the logarithm of negativity .55J to 
accentuate small variations from E/j = 0. 



B. NPT bound for {N/2, N/2} bipartite split 

The values of N and a for which the minimum eigen- 
value of p^ vanishes give a bound on entangled p C H- In 
what follows, we show that it is possible to find an ana- 
lytic expression for the minimum eigenvalue in terms of 
N and a. 

Numerical calculations reveal that at fixed N , the min- 
imum eigenvalue A m i n appears to correspond to one of 
two eigenvectors depending on whether a is above or be- 
low a transition value. Based on this evidence, we conjec- 
ture that the eigenvector corresponding to the minimum 
eigenvalue is 



\v-), a < a tr 
\v + ), a > a tr 



(33) 



with at r being the polarization where the transition in 
eigenvectors occurs. The eigenvectors are given by 



\v-) 



1 



i) 



1 



| 2 Ar-r 
12* - 2 N / 2 ) 



(34) 
(35) 



in the computational basis. The state labels are un- 
derstood in binary notation. We numerically verified 
Eqs. up to N = 10. 

We now obtain an analytical formula for the hypoth- 
esized Amin- Using the relation pj^\v±) = X±\v±), we 
find 



2 N~1 _ 1 



X- = 

-(2^-1 

A + = 6 w/2 -i 

- 2 N / 2 _ 



rcH 



2 N-1 _ 1 



o Ta 



2 N ' 2 - 1 



o Ta 



2^ _ 2 N / 2 



(36) 



(37) 



The minimum eigenvalue is A m i n = min(A_, A + ). 

The matrix elements in Eqs. (|36[1 and 1)37(1 are cal- 
culated in the following manner. Observe that a given 
matrix element {i\p^.\j) = (k\p c n\£) where \k){£\ — 
(I*)0I) Ta - Since U c n is sparse and p tri is diagonal in 
the computational basis, (k\p c n\£) may be expressed as 
a simple sum of diagonal entries in p t h- Formulas for the 
relevant matrix elements of are given in Table [I] 

Substituting the results from Table Q] into Eqs. I|36|)- 
l(37|l and applying Eq. H17|) . we obtain 



e Q (e cosha — sinha) 



e "(cosh a — e sinha) 



(38) 
(39) 
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{N/2.N/2} split 
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iog 10 (a _1 ) 



FIG. 3: (color) Bounds on nonseparable CNOT-Hadamard- 
transformed thermal states under {l,iV — 1} bipartite split 
(Eq. and {N/2, N/2} bipartite split (Eq. 0TJ). The 

light red region marks parameter space where p C H is nonsep- 
arable under the {N/2, N/2} split. The union of the light red 
and light blue regions shows where p ct i is nonseparable un- 
der the {1, N — 1} split. Effective pure states are entanglable 
in parameter space left of the Braunstein et al. lower bound 
(Eq. 1121 ') and are nonentangiable in the parameter space be- 
neath the Gurvits-Barnum bound (Eq. mai l. 



The transition value atr is determined by the condition 
A + — A_ =0 and can be calculated by solving the relation 



e Nat * tanh 



1 



(40) 



The expressions for A± show that Eq. I|40|l is also a 
constraint for A m i n = 0. Therefore, the bound on non- 
separable p C H can be specified by 



,JVct 



tanh a > 1 



Discussion 



(41) 



We compare the NPT bounds on the entanglability 
of thermal states with comparable bounds on effective 
pure states in N-a parameter space. Fig. [3] shows the 
bounds on nonseparable p cr i for the {1,N — 1} and 
{N/2, N/2} split in N-a space. In both cases, the ther- 
mal state bounds include more parameter space than the 
Braunstein et al. lower bound, although the {1,N — 1} 
split gives the tighter constraint. Neither thermal state 
bound succeeds in capturing parameter space where p e s 
is nonentangiable. 

Numerically solving Eq. I|31|l . we find that for polar- 
ization a = 4 x 10~ 5 , we require at least 126,584 qubits 
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Matrix element Equivalent expression 

(2"- 1 - ljp^- 1 - l) i((2 JV - 1 - ljp^- 1 - l) + (2* - l\p th \2 N - l)) 

^-i-ll^^- 1 ) i((2 JV - 1 -2 JV / 2 | /5th |2 JV - 1 -2^) -(2 JV -2 JV / 2 |p th |2 JV -2 JV / 2 )) 

(2^ - llp^f 12^ - 1) ^((a^ 2 - llpthh^ 3 - 1) + (2^ + - lLhb^ 2 + 2 JV ~ 1 - 1)) 

(2^ 2 - l\p%\2 N - 2^ 2 ) |((o|pth|o) + (2 N ~ 1 \p th \2 N - 1 )) 



TABLE I: Formulas for several matrix elements of p c £ given in terms of the diagonal elements of p t h- These expressions are 
needed to derive an analytical formula for the minimum eigenvalue of p^ under the {N/2, N/2} bipartite split. 



to achieve entanglement for ?7 c h transformed thermal 
states. This number is clearly impractical, but it dra- 
matically improves upon the Schulman-Vazirani value of 
6.3 x 10 s qubits from Sec.lTTTl 



VI. BOUNDS ON THE SEPARABILITY AND 
DISTILL ABILITY OF BELL-TRANSFORMED 
THERMAL STATES 

We now apply a more general procedure to derive 
bounds on the separability and distillability of Bell- 
transformed thermal states. Unlike the methods of 
Sec. [V] this procedure allows us to find bounds on the 
entanglement of thermal states transformed by any Bell 
unitary under any bipartite split. The main idea is to 
take a state p' with known entanglement properties and 
find a random local unitary operation that transforms 
UbPthUl 1— > p' . Because such an operation cannot in- 
crease the entanglement of the initial state, the entangle- 
ment of pf gives a lower bound on the entanglement of 
Bell-transformed thermal states. 

A candidate state for p' is provided by a family of 
mixed Bell states Pn whose entanglement has been clas- 
sified by Diir and Cirac [24|. We first review their for- 
malism and then apply it to Bell-transformed thermal 
states to find bounds on their separability and distilla- 
bility. These results give tighter bounds than the direct 
calculations of Sec. ]V\ 



A. Diir-Cirac classification of entanglement in 
special mixed Bell states 

1. Formalism 

Diir and Cirac consider a special family of mixed Bell 
states parameterized by the number of qubits N . They 



are denned by 

PN = A+|*+)<*+|+A -|* -)<*ol 

2 («-i>_i 

3=1 

where the generalized Bell states are given by 

|*±> = ^(|0>|j>±|l)|j>) (43) 

in the computational basis. The state of the first qubit 
and the state of the latter N—l qubits are explicitly sep- 
arated |Ei|- Here \j) = \2 N ~ 1 - j - 1) is the bit-flipped 
version of \j), and 1 < j < 2 W ~ 1 — 1. The parameters 
A q are chosen such that A^ — Aq~ > 0. 

We specify a bipartite split with a nonnegative integer 
k, which when expressed in binary, labels the qubits that 
are in party A with Is and the qubits in party B with 
0s. For example, if k = 010110, then the second, fourth, 
and fifth qubits arc in party A, and the first, third, and 
sixth qubits are in party B. With no loss in generality, 
we require the first qubit to always be in party B, i.e. the 
most significant bit is always 0. Therefore, the possible 
bipartite splits lie in the range 1 < k < 2 1 — 1. 

The family of states pn has particularly elegant con- 
ditions on its partial transposition. A straightforward 
calculation shows that pn under bipartite split k has the 
properties: 

A < 2X k ^ p N has PPT (44) 
A > 2X k p N has NPT (45) 

where we have defined A = Aq — Aq~ . 

More importantly, we can characterize the full separa- 
bility and distillability of pn- Diir and Cirac proved the 
following: 

1. Consider all possible bipartite splits of an N qubit 
system. If and only if each of these splits has PPT, 
then pat is fully separable, i.e. separable under any 
partition of the system. 
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2. Consider all possible bipartite splits of an N qubit 
system. If and only if each of these splits has NPT, 
then pat is fully distillable, i.e. a maximally entan- 
gled pair can be distilled from any two particles in 
the system. 

Thus the PPT and NPT conditions give bounds on the 
separability and distillability of pn- 

A < 2 minjfcj (Afc) pn fully separable (46) 
A > 2 max{ fe j (Afc) => pn fully distillable (47) 

where the minimum and maximum are taken over all 
possible bipartite splits k. 

The bound on full distillability is significant because 
it specifies the regions in parameter space where useful 
entanglement can be obtained in the form of maximally 
entangled pairs. 



2. Bound on the entanglability of effective pure states 

As Diir and Cirac discuss [24|, the fully distillable cri- 
terion gives a simple bound on the entanglability of ef- 
fective pure states that is tighter than the Braunstein et 
al. lower bound in Eq. 0. 

Suppose we apply a unitary to the effective pure state 
that transforms |0) to the maximally entangled Bell state 



I* 



/ , yielding a new state 
1-e, 



Pcff = 



(48) 



that fits the form oi p^. It is easy to see that Xq = 
(l-e)/d + e and A = Xj = (l-e)/d. Therefore, A = e. 

Eq. (gTJ) shows that full distillability of the trans- 
formed p c ff requires 



e > 



(49) 



1 + 2*- 1 ' 

Substituting Eq. i|ll|) . we obtain a bound on the en- 
tanglability of effective pure states in NMR parameter 
space: 



a > 



1 



In 



2 + 2 



iY 



1/N 



(50) 



B. Application to Bell-transformed thermal states 

1. Random, local unitary operation mapping 
Bell-transformed thermal states to pjv 

The thermal state is diagonal in the computational ba- 
sis, so any Bell-transformed thermal state has the form 

Pb = A+|*t)<*+| + A^|*7)<*7| (51) 

2 («-i>_i 

+ e (wwi+ww 



The state is diagonal in the generalized Bell basis 
{l^j)}> Du t i s n °t °f the form p^ because generally 
A+^A-. 

Suppose we can find a random local unitary operation 
£g such that £i(pb) has form p^. The action of any ran- 
dom local unitary operation on an arbitrary p can always 
be expressed as a probabilistic mixture of unitaries Uc 



£i(p) = Y.p* u *P u > 



(52) 



where pi is the probability of applying the unitary Ui. 
All the unitaries Ui are local, meaning that they act 
on each qubit separately. Therefore the entanglement 
of each state UipU\ cannot be larger than the entangle- 
ment of the original state p. Moreover, the quantum op- 
eration £i(p) can be interpreted as a convex combination 
of transformed thermal states. Thus the entanglement of 
£e{p) must be less than or equal to the entanglement of 
p itself; a random mixture of local unitaries cannot in- 
crease entanglement. It follows that the entanglement of 
£i(pb) must bound the entanglement of ph- 

Now we describe a particular random local unitary 
operation that performs the desired transformation. 
Consider the following procedure: 

1. Start with the Bell-transformed thermal state ph- 

2. Apply mixing operation R = Ri where Ri 
multiplies the zth qubit by a random phase 4>i if the 
qubit has state |0) and Ri does nothing otherwise. 
We assume that <pi is uniformly distributed over 
[— 7T, 7r] subject to the constraint £\ <fn — 2n. This 
requirement is chosen so that R\^q) = l^o)- 

3. Average Rph over all <f>i. 

Let us verify that £§ produces a state of form p^. 
First, we calculate the effect of R on each generalized 
Bell state: 



\xf) 



^(^|0)|j)±e-^|l)|j)) ,5+ 



where 8j, the phase of |0)|j) due to the action of R, is 
given by 



W(m)=i} 



(54) 



The function / is unity only if the «th binary digit of j 
is one. Note that |l)|j) acquires a phase of —0j because 
it is the bit-flipped version of |0)|j). 

Next, we find the state that results when R acts on the 
entire Bell-transformed thermal state: 

2 (N-i)_ 1 

R Pb tf= J2 (*t\xt){4\+*J\xJ)<Xj\) (55) 
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where 

\xf){xt\ 



1 - 
2 

± 



(|0)|j)(0|(j| + |l)|j)(l|(j|) (56) 
(e 2 ^|0)|j)<l|(j| + e ^|l)|j)(0|<j| 



When we average Eq. Q55[) over fa, the last two terms 
in Eq. ljSl)jl vanish for all j ^ 0. When j — 0, E^, has no 
effect as 6 is constrained to be 2ir. The final state after 
application of E^ is consequently 

E^(p h ) = A+|*+)(*+| + Ao|* -)(*ol (57) 

2<~- 1 >-i A + + A - 

which is indeed of form pj^. Matching this expression 
with Eq. (@2J>j we obtain 



A, 



(58) 



Inserting the above into Eqs. I|44|) and l|45(l . we estab- 
lish PPT and NPT conditions on p c n under bipartite split 
k: 

tanha < e ' 2w[k)a => p cH has PPT (64) 
tanha > e r 2w{ ^ a =>. PcH h as NPT . (65) 

From the above conditions, the entanglement of /?b is 
constrained by the bounds 

tanha < e" 2(JV - 1)a ^ p b fully separable (66) 
tanh a > e~ 2a p h fully distillable . (67) 

The separable bound corresponds to the {1,N—1} split, 
and the distillable bound corresponds to the {N — 1,1} 
split. 



Separability and distillability of U c HUf an transformed 
thermal states 



The effect of is to average the eigenvalues A^ . 

We emphasize that is not the only possible opera- 
tion that produces the needed transformation. Random 
local unitary operations generally decrease the entangle- 
ment of the system. We desire a procedure that preserves 
as much entanglement as possible. In the next subsec- 
tion, we will make a few remarks about the effectiveness 
of £,f, in this regard. 



2. Separability and distillability of U c h transformed thermal 

states 



We now calculate separability and distillability bounds 
on Bell-transformed thermal states. First, we exam- 
ine the case where the Bell unitary is chosen to be the 
CNOT-Hadamard transformation U c n. 

Fig-fflshows that [/ c h maps computational basis states 
to Bell states in the following manner: 



I* 



)N-1 
iN-l 



i-2(* 



N 



Inserting Eq. (|17fl . we find 

At = (*>cH|*t) 



A; 



,[N-2w(j)]a 



,[N-2w(j)-2]a 



(59) 



(60) 
(61) 



Then we apply the random local unitary operation E^ 
and find that E^Pca), when expressed in the form of pn, 
has parameters 



1 



e Na 1 _ e -2« 



(62) 



A fc = ±e^ N - 2w ^ a {l + e' 2a ). (63) 



We turn to the case where the Bell unitary is chosen to 
be U c nU[ an . This unitary yields tighter bounds compared 
to £/ch, as we now show. 

The tightness of a bound derived from Diir-Cirac for- 
malism increases with the size of A. In the case where we 
transform p t h with U c u, we effectively set Xq — e Na /Z 
and Aq = e^ N ~ 2 ' a /Z, giving a gap A that scales as 2~ N . 
However, a larger gap may be achieved. 

We can shuffle the unitary mapping between compu- 
tational basis states and Bell states to obtain the values 
of A q and A.,- that maximize A. The largest possible gap 



occurs when we select An 



„±Na 



jZ, which yields 



A = 



2 sinh(TVa) 



(68) 



This choice scales as N2 N . 

The new gap can be realized with mapping 



< i < 2 



N-l 



1 



< i < 2 1 



(69) 



This mapping is exactly produced by the Bell unitary 
UcnUten, with the permutation matrix [/f an shuffling the 
old mapping in Eq. |1I1> to give Eq. 

We proceed to find bounds on the separability and dis- 
tillability of pb- Doing the same calculation as before, we 
have 



A t 



1 



cosh 



(jST-2w(kfja . (70) 



Comparing Eqs. I|68|l and (|70|l . we establish the follow- 
ing conditions on p c { under bipartite split k: 



tanh(fca) < e -W-v>W]'* 
tanh(fca) > e -^ N ~ w {k)]a 



Pd has PPT (71) 
p cf has NPT . (72) 
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FIG. 4: (color) Bounds on the separability and distillability of Bell-transformed thermal states as derived from Dur-Cirac 
formalism, (a) 2 < N < 20 and (b) 2 < N < 500: p cH fully separable (Eq. JSHJ), p C H fully distillable (Eq. (JSTJ), Pcf fully 
separable (Eq. (17310 . and p c f fully distillable (Eq. 11741 '). Transformed thermal states p C H and p c { are separable/distillable in 
the parameter space to the right /left of the appropriate separable/distillable bound. Effective pure states are nonentanglable 
in the parameter space beneath the Gurvits-Barnum bound (Eq. 1131 ) and are entanglable in the parameter space left of the 
Dur-Cirac bound (Eq. (IBUl). 



Minimizing and maximizing the above conditions over 
all possible bipartite splits, we find that py, is fully sepa- 
rable if and only if 



sinh(iVa) < 1 
and that pb is fully distillable if and only if 
tanh [(N 



-2a 



l)a\ > e 
tanh a > . 



(73) 



(74a) 
(74b) 



Eq. O is der ived from the {N/2, N/2} split whereas 
Eqs. I|74a|) and (|74bjl correspond to the {1,N - 1} and 
{N — 1, 1} splits respectively. These two splits give the 
same bound because in Eq. (|70|l is symmetric under 
the replacement oi w(k) by N — w{k). 



C. Discussion 

We first evaluate the effectiveness of our random lo- 
cal unitary operation E^. Observe that the {1,N — 1} 
and {N/2, N/2} bipartite splits correspond to Hamming 
weights w(k) of N — 1 and N/2 respectively. If these 
weights are inserted into Eq. (|65|l . we recover the NPT 
bounds calculated in Eqs. I|31f) an d (|41|l . We also nu- 
merically computed NPT bounds (up to N — 10) for 
the bipartite splits corresponding to the bounds on fully 



separable and distillable pb in Eqs. (|rlfj|) . l|rj7|) . l(75|) and 

(|74|l . The excellent agreement in all cases implies that 
our procedure £ § does not erase any of the information 
we obtain from the direct NPT test. 

We now examine the behavior of the Dur-Cirac de- 
rived bounds on fully separable and distillable Bell- 
transformed thermal states in NMR parameter space. 
These results are shown in Fig. 21 As a comparison, we 
have also plotted bounds on the entanglability (Diir and 
Cirac) and nonentanglability (Gurvits and Barnum) of 
effective pure states. Note that the bound on entanglable 
effective pure states also implies full distillability of trans- 
formed Poff- 

The U c nUf an unitary produces tighter bounds than 
C/ c h, as we expect. Curiously, the bound on fully dis- 
tillable /9 c h exactly coincides with the Braunstein et al. 
bound on the nonentanglability of p e s (see Fig- , and 
the bound for fully separable p C H exactly matches the 
bound on fully distillable p c {. 

We are interested in knowing whether our thermal 
state bounds improve upon the known constraints on ef- 
fective pure states. First, we compare the bound on fully 
distillable p c { against the corresponding bound for effec- 
tive pure states. The thermal state bound encompasses 
more parameter space at small N, but the effective pure 
state bound does better for N > 9. Yet thermal states 
should be at least as entanglable as effective pure states. 
The apparent contradiction indicates that while U c nUfa n 
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is the optimal Bell unitary, it is not the optimal unitary 
in general for entangling thermal states. Second, we ob- 
serve that the Dur-Cirac derived bounds on separable pb 
are much looser than the Gurvits-Barnum bound on the 
nonentanglability of p c ff. This result again suggests that 
the family of unitaries Ub is not optimal. 

We do succeed in finding entangled thermal states out- 
side the Dur-Cirac bound on entanglable effective pure 
states. Thermal states outside the parameter space of 
fully separable pb have at least one NPT bipartite split, 
implying that at least one maximally entangled pair can 
be distilled. The bound on fully separable p c f (Eq. J73J)) 
thus gives the tightest bound on entanglable thermal 
states, scaling as N ~ 1/a. It also captures new pa- 
rameter space outside of the Dur-Cirac bound on en- 
tanglable p c g. Assuming the largest realistic initial po- 
larization a = 4 x 1(T 5 , Eq. O requires N < 22,035 
for entanglable thermal states. In comparison, the Dur- 
Cirac bound in Ea.1501 gives a constraint N < 50, 000 for 
entanglable effective pure states in the same parameter 
space. 



VII. MAJORIZATION APPROACH TO 
DETERMINING ENTANGL ABILITY 

We have seen that the Dur-Cirac formalism is a power- 
ful tool for obtaining bounds on the entanglement of Bell- 
transformed thermal states. The key is to find a random 
local unitary operation that maps the Bell-transformed 
thermal state to a state of form pjy. In principle, we 
can always try to construct the needed quantum oper- 
ation, but this method is time-consuming and wasteful. 
Knowing such a procedure exists is enough. 

Motivated by the connection between entanglement 
and the mathematical theory of majorization |47|. we 
discuss a different approach that uses majorization to 
address this problem. We first define majorization 
and describe two theorems from majorization theory: 
Uhlmann's theorem and the majorization condition on 
von Neumann entropy. Then we explain how these the- 
orems may be used to potentially derive tighter bounds 
on entanglable thermal states. 



A. Formalism 

Majorization is a mathematical relation that states 
whether one probability distribution is more disordered 
than another. 

Suppose we have two probability distributions that 
can be described by real vectors of length d; r = 
(n, rs, ■ ■ ■ r<j) and s = (si, S2, . . . , Sd)- Let and s' be 
the vectors r and s with their components arranged in 
non-decreasing order, that is r{ < r 2 < ... < r d and vice 
versa. We say "s majorizes r" or "r is majorized by s" if 



and only if 

X>J>£4 (75) 

3=1 3=1 

for k = 1, 2, . . . , d and with the equality holding for k = d. 
This relation is denoted by r -< s. When it holds, we may 
think of r as being more "disordered" than s [4|| ■ 

In the context of quantum mechanics, we take r and s 
to be vectors containing the eigenvalues of density matri- 
ces. Given two density matrices p and cr, we have p -< a if 
and only if A(p) -< A(er) where A(p) is a vector containing 
the eigenvalues of p. This interpretation makes sense be- 
cause the eigenvalues define a probability distribution on 
an orthogonal set of states. Thus, majorization for two 
density matrices is equivalent to comparing the disorder 
of their corresponding probability distributions. 

There are two results from majorization theory that 
are useful for our purposes: Uhlmann's theorem |49i IHoL 
l5lll52j | and the majorization condition on von Neumann 
entropy [Igj . 

Uhlmann's theorem states that given two Hermitian 
matrices A and B, there exists a random unitary opera- 
tion £ such that A — £{B) if and only if B majorizes A. 
Mathematically, we write 

A -< B 3 £ s.t. A = £ (B) = ^piUiBUj . (76) 

i 

Note that unlike the operation £g of Sec. lVIl the unitaries 
Ui here need not be local. 

The majorization condition on von Neumann entropy 
S is the following. Given two density matrices p and a, 
if a majorizes p, the von Neumann entropy of p must be 
at least as large as the von Neumann entropy of a, i.e. 

p ~< a =► S(p) > S{a) . (77) 

Notice that the majorization relation is merely necessary 
and not sufficient. 



B. Application to thermal states 

Now we explain how the formalism we just described 
may be used to derive bounds on the entanglability of 
thermal states. 

Suppose pth majorizes p', a state of form pjv. Then 
Uhlmann's theorem says that p' = X^i^AhC^/- The 
same convexity argument from Sec. IIIII establishes that 
the entanglement of p' gives a lower bound on the maxi- 
mum entanglement attainable from thermal states. 

We have converted the problem of calculating the en- 
tanglement of U ptbU^ to finding a density matrix p' that 
is majorized by pth- Supposing that p' is well-defined in 
terms of N and a, we can derive analytical bounds on 
entanglable thermal states. 

Unfortunately, this majorization approach is not nec- 
essarily simpler. We need to guess the right form for p', 
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a state with 2 + 1 parameters. The only constraint is 
that p' be entangled with respect to at least one bipartite 
split. Then we must compare as many as 2 N eigenvalues 
to see if the majorization actually holds. 

However, we can use Eq. (|77|l as a first cut criterion 
before checking that the desired majorization relation 
holds. When S(p t h) > S(p'), it is impossible for p t h to 
majorize p' . The von Neumann entropy of the thermal 
state has the simple formula 

%) = -JV[plog 2 (p) + (l-p)log 2 (l -p)} (78) 

where p = (1 + tanha)/2. So, if there is a closed form 
expression for S(p'), it is easy to see if Eq. (|77|l holds. 
This method may save time in searching for the best 
construction of p' . 

We propose two different strategies. To simplify mat- 
ters, we hold a constant. The first idea is to guess 
an analytical form for p', such that at least one bipar- 
tite split has NPT. We can derive a formula for S(p') 
and use the majorization condition on von Neumann en- 
tropy to find the minimum number of qubits needed for 
S(pth) > S(p'). We then set N to this minimum value 
and increment iV until the majorization relation p' -< pth 
holds. Alternatively, we can attempt to numerically op- 
timize p' to be majorized by /?th at fixed N, increasing 
N if the optimization fails. 

While we have not yet found a good construction 
for p', we believe this majorization-based approach may 
be fruitful in future studies of NMR thermal state en- 
tanglability. 

VIII. CONCLUSION 

In summary, we have obtained a bound on the en- 
tanglability of thermal states in liquid-state NMR by 



classifying the entanglement of Bell-transformed thermal 
states. These results are confirmed by more direct calcu- 
lations of NPT bounds for two specific bipartite splits. 



This paper sought to investigate two related issues: 1) 
the nature of mixed state entanglement in quantum com- 
putation and 2) how to most efficiently achieve entangle- 
ment in a liquid-state NMR system. We have made some 
progress on both accounts. Our tightest lower bound on 
thermal states scales as N ~ 1/a and provides us with 
a sense of the resource tradeoffs that must be made to 
achieve entanglement. Moreover, this result explicitly 
shows that thermal states are more easily entangled than 
effective pure states. Our bound identifies entanglable 
thermal states in N-a space outside the best bound on 
entanglable effective pure states. Assuming the largest 
reasonable initial polarization (a = 4 x 10~ 5 ), we find 
that a thermal state must have N < 22,305 to be en- 
tanglable, a limit that improves upon the Dur-Cirac con- 
straint of N < 50, 000 for effective pure states. 



Our methods exploit the symmetry of Bell unitaries, 
but there is considerable room for improvement, espe- 
cially if these studies of liquid-state NMR entanglement 
are to motivate practical experiments. We have sug- 
gested a majorization-based approach that uses a state p' 
of form pjv to bound the entanglability of thermal state 
without having to choose a specific unitary. The missing 
piece is an efficient construction of p' — a task we leave 
for future work. 
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